Temperature dependence of Vortex Charges in High Temperature Superconductors 
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Using a model Hamiltonian with d-wave superconductivity and competing antiferromagnetic (AF) 
interactions, the temperature (T) dependence of the vortex charge in high T c superconductors 
is investigated by numerically solving the Bogoliubov-de Gennes equations. The strength of the 
induced AF order inside the vortex core is T dependent. The vortex charge could be negative when 
the AF order with sufficient strength is present at low temperatures. At higher temperatures, the 
AF order may be completely suppressed and the vortex charge becomes positive. A first order 
like transition in the T dependent vortex charge is seen near the critical temperature Taf- For 
underdoped sample, the spatial profiles of the induced spin-density wave and charge- density wave 
orders could have stripe like structures at T < T s , and change to two-dimensional isotropic ones at 
T > T s . As a result, a vortex charge discontinuity occurs at T s . 

PACS numbers: 74.20.-z, 74.25Jb 



For the past many years there have been intensive 
studies of the vortex physics in high temperature super- 
conductors (HTS). It is well established now that the 
d-wave superconductivity (DSC) with a competing an- 
tiferromagnetism (AF) order plays an important role in 
determining the vortex structure of HTS. Many theo- 
retical studies [J H H 0, H, H 0, H H [HI have shown 
that the AF order may appear and coexist with the un- 
derlying vortices. Experimental facts including neutron 
scattering [11] , moun spin rotation measurement Il2l and 
nuclear magnetic resonance (NMR) experiments |lA 
have provided a strong support for the existence of AF 
order inside the vortex core in appropriately doped HTS. 
On the other hand, the electrical charge associated with 
the vortex in superconductors has also been paid consid- 
erable attention both theoretically [H Hill E E 
and experimentally |2^,|23. Based on the BCS theory, 
Blatter et al. pointed out that for s-wave supercon- 
ductor the vortex charge is proportional to the slope of 
the density of states at the Fermi level. However, the 
NMR and nuclear quadrupole resonance (NQR) mea- 
surements on YBa2Cu307 and YBa2Cu408 |2l| seemed 
to obtain results for the vortex charge, contradictory to 
that predicted from the existing BCS theory with respect 
to both sign and order of magnitude. In addition, if the 
impact of vortex charge on the mixed state Hall signal is 
considered ^| , the sign estimated from Hall effect exper- 
iments |22| for various HTS materials disagrees with the 
prediction of BCS-type theory. In view of these points, 
together with the fact that the strong electron correlation 
with the (i-wave superconducting pairing was not consid- 
ered for the vortex charge, we carried out the previous 
studies [2(J at zero temperature and found that the vor- 
tex charge in HTS is strongly influenced by the presence 
of the induced AF order in the vortex core. The charge 
carried by a vortex is always positive for a pure d-wave 
superconductor, and it becomes negative when there is 
sufficient strength of AF order in the vortex core poj . 



Since most relevant experimental phenomena have been 
observed at finite T not close to zero, including the sign 
change of the mixed state Hall effect it is necessary 
and valuable to study theoretically the vortex charges 
at finite T, which so far has not been addressed in the 
literature. 

In this article we shall investigate the vortex charge 
as a function of T for both optimally doped and under- 
doped HTS in detail. Our calculation will be based on 
the model Hamiltonian on a square lattice with a nearest 
neighboring attractive interaction Vdsc describing the 
DSC and an onsite Coulomb repulsion U representing the 
competing AF order. We shall adopt the well-developed 
numerical scheme to study the vortex structure. Our 
numerical results at finite T confirm that the sign of the 
vortex charge is still determined by the AF order in the 
vortex core, being consistent with our previous work poj . 
The vortex charge could be transformed from negative 
to positive values as T is increased to a critical value 
Taf such that the AF order is completely suppressed at 
T > Taf-, and the transition seems to be first-order-like. 
For underdoped sample we show that the symmetry of 
the induced spin-density wave (SDW) and the charge- 
density wave (CDW) may change from stripe like to two 
dimensional like at a critical temperature T s . A discon- 
tinuity of the T-dependent vortex charge also appears at 
T a . 

The effective mean-field Hamiltonian describing both 
DSC and SDW orders on a two-dimensional lattice can 
be represented as, 

H = -^UjctrCja+^iUriia - [i)c\ G C i(7 

where c\ a is the electron creation operator, \i is the chem- 
ical potential. In the presence of magnetic field B per- 
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pendicular to the plane, the hopping integral can be ex- 
pressed as = £oexp[z^- j* % A(r) • dr] for the nearest 
neighboring sites A Landau gauge A = (—By, 0, 0) 

is chosen. The two possible SDW and DSC orders in 
cuprates are defined as Af DW = U{c\^c^ — c[^c%\) and 
Aij = V DS c(c^Cji ~ Q|Cj|)/2, where U and V DS c rep- 
resent respectively the interaction strengths for two or- 
ders. The mean-field Hamiltonian (1) can be diagonal- 
ized by solving the resulting Bogoliubov-de Gennes equa- 
tions self-consistently 

e( A*f -nij{lf;)= E -( u ^;) - ( 2 ) 

where the single particle Hamiltonian Hij^ = — Uj + 
(Un i& - n)5ij, and n iT = J2n\ u i^\ 2 f( E n), ri u = 

En W 1 ~ f(E n ))> \j = ^En«T^* + 

v*^ ) tanh ^ 2 ^ T ^ , with /(£") as the Fermi distribu- 
tion function. The DSC order parameter at the ith site 
is Af = (Af +e ^ + Af_ ex>i - AP i+e> - A£_J/4 where 

= A^exp[z^ J^ r " +r ^/ 2 A(r) • dr] and e X)2y denotes 
the unit vector along (x, direction. In our calcula- 
tion, the related parameters are chosen as: t = a = 1, 
Vdsc = 1-2, the linear dimension of the unit cell of the 
vortex lattice is TV^ x iv"^ = 40 x 20. This choice corre- 
sponds the magnetic field B c± 37T. 

In the absence of the magnetic field and for optimally 
doped sample (S = 0.15), the parameters are chosen in 
such a way that only DSC (with T c ~ 0.164) prevails in 
the system. Our numerical results with U = 2.2 indeed 
show that the AF order is induced at low temperature 
and is suppressed when T is larger. The spatial profiles of 
the superconductivity order parameter, staggered mag- 
netization, and the charge density near the vortex core 
are plotted at two different temperatures in Fig.l. Here 
the staggered magnetization of the induced AF or SDW 
order is defined as M] = (-l)*Af DW /U. The AF order 
is generated in the region where the DSC order parame- 
ter is partially suppressed. The low temperature results 
at T = 0.02 are shown in Figs. 1(a) to 1(c) where the AF 
order is nucleated and spreads out from the core center, 
and the vortex charge is negative. The induced AF order 
behaves like a two dimensional SDW with the same wave- 
length in the x and y directions. The higher temperature 
results at T = 0.05 are shown in Figs. 1(d) to 1(f) where 
the AF order is completely suppressed and the vortex 
charge becomes positive. The appearance of the SDW 
order pinned by the vortex lattice at low T strongly en- 
hances the net electron density (or depletion of the hole 
density) at the vortex core as shown in Fig. 1(c). The 
present result is consistent with that of Ref. 20 where 
the doping and U-value dependences of vortex charges 
are examined. Both works reach the same conclusion: 
the vortex charge is negative when a sufficient strength 
of SDW order is induced , while it turns to positive with 
the suppression of SDW order due to increasing T or 
doping level, or decreasing the [/-value. 



Fig. 2(a) plots the phase diagram of T versus U 
for positively (hole-rich) and negatively (electron-rich) 
charged vortices. It is obvious that the AF vortex core 
can easily be induced in the low temperature regime or 
with stronger AF interaction U while the pure DSC core 
tends to exist in the high T or smaller U . The electron 
density inside the core is higher than the average density 
in the low temperature region while the electron density 
becomes lower than the average in the high temperature 
region. There exists a clear boundary between these two 
phases. To estimate the vortex charge, we first deter- 
mine the vortex size by examining the spatial profile of 
DSC order parameter. Next we sum over the net elec- 
tron density inside the vortex core. In Fig. 2(b), we 
depict the chemical potential versus T at fixed doping 
level (5 = 0.15 with U = 2.4 and U = 2.2). The result 
(solid line) for U = 2.4 exhibits a first-order like transi- 
tion at Taf ~ 0.077, indicating the existence of induced 
AF order below this critical temperature. For U = 2.2, 
the result (dashed line) seems to show a very weak dis- 
continuity which can hardly be seen in Fig. 2(b). The 
above finding may be closely related to the experimen- 
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FIG. 1: Spatial variations of the DSC order parameter 
[(a) and (d)], staggered magnetization Mf [(b) and (e)], and 
net electron density m [(c) and (f)] in a 20 x 20 lattice. The 
left panels [(a), (b), and (c)] and the right panels [(d), (e), 
and (f)] are for T = 0.02 and T = 0.05, respectively. The 
averaged electron density is fixed at n — 0.85 and U — 2.2. 
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FIG. 2: (a) Phase diagram of interaction strength U versus 
temperature T for positive and negative charged vortex at 
optimal doping 5 = 0.15. (b) temperature dependence of 
chemical potential for U = 2.4 and U = 2.2. 
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FIG. 3: The temperature dependence of the number of vortex 
charge Q v /e for U = 2. 4 (solid line), U = 2.2 (dotted line) and 
U = (dashed line). The averaged electron density is fixed at 
ft = 0.85. 



tal results for underdoped TlBa2Cu408 where the 
vortex core with induced AF order below Taf was re- 
ported. Also from Fig. 2(b), a discontinuity in the slope 
of the T-dependent chemical potential at Taf for both 
cases can be clearly observed. We expect that the local 
density of states (LDOS) of the vortex core should have 
double peaks 0, H Efl 123 13 near zero °i as for T < Taf', 
while for T > Taf only a single broad zero-bias peak can 
appear in the LDOS as one of the essential character- 
istics of DSC [22|. It would be interesting to measure 
this temperature evolution in the LDOS by future STM 
experiments. Also interestingly, even though the origin 
of Hall sign anomaly in HTS is still debatable 24], the 
vortex charge could make an additional contribution to 
the sign change in the mixed state Hall conductivity . 
As a result, we predict that a sign reversal in the Hall 
signal may occur at T not too far below Taf- 

In Fig. 3, we plot the vortex charge number Q v /e as 
a function of T for three typical [/-values, here e < 
is the charge of an electron. The T dependence of the 
induced staggered magnetization M£ at the vortex core 



center has also been examined in previous studies [53, 
where the induced AF order is shown to disappear for 
T > Taf- From Fig. 3 and for U = 2.4, one can clearly 
see an abrupt jump in the number of vortex charge Q v /e 
as T varies around Taf ~ 0.077 which is consistent with 
the critical temperature as shown in Fig. 2(b). This 
positive to negative vortex charge transition seems to be 
first-order like. The magnitude of the discontinuity re- 
duces to one third when U = 2.2 at Taf ~ 0.04. Here 
the discontinuity occurs at the positive to positive vortex 
charge transition, not the positive to negative transition 
as studied in the case of U = 2.4. This is because slightly 
below Taf, the induced AF order is too weak to make the 
vortex charge negative. For a pure DSC case with [7 = 0, 
the vortex charge Q v is always positive and its magnitude 
first increases slightly with T and then decreases to zero 
as T approaches to T c . The critical temperature Taf is 
U- value dependent or the sample-dependent. The larger 
U corresponds to larger Taf- These results confirm that 
the vortex charge is strongly influenced by two compet- 
ing effects: the suppression of the DSC order at the core 
center leading to the depletion of the electrons, and the 
induction of the AF order which favors the accumulation 
of electrons. The condition for the negative vortex charge 
to appear depends solely on whether there is a sufficient 
AF order inside the vortex core. We emphasize that our 
results are robust and indifference to band parameters, 
and should give a qualitative description on the vortex 
physics in HTS. 

In the following we are going to investigate the tem- 
perature dependence of the vortex charge for underdoped 
HTS (6 = 0.12). For U = 2.4, the spatial profiles of the 
staggered magnetization and electron density near the 
vortex core at three different temperatures are plotted 
in Fig. 4. As discussed in our previous work 7], the 
parameters used to obtain Fig. 4 yield stripe like struc- 
tures in the DSC, SDW and CDW order parameters at 
very low temperature. In the present work, the DSC 
and SDW orders coexist in underdoped sample below T c , 
which implies Taf > T c . So the question arises: will 
these stripe like structures persist with increasing tem- 
perature? From our calculations, the spatial profiles of 
these order parameters may go through three different 
symmetries as T increases. At T = 0, the spatial distri- 
bution of SDW (M s ) and CDW (ni) orders exhibit quasi- 
one-dimensional stripe like behavior (see Figs. 4(a) and 
4(b)) . When T is raised to T = 0.11 above a characteris- 
tic temperature T s (~ 0.1), the symmetry of the patterns 
change to isotropic and two-dimensional (see Figs. 4(c) 
and 4(d)). Above T c (~ 0.16) and at T = 0.17, the resid- 
ual AF and electron density become uniform as shown 
in Figs. 4(e) and 4(f). The AF order vanishes when 
T > T AF - 

Finally we display the temperature evolution of the 
vortex charge number Q v /e for the underdoped case 
S = 0.12 with U = 2.4 and U = 2.2 in Fig. 5. Comparing 
with Fig. 3, the magnitudes of vortex charges for both 
cases in Fig. 5 at low temperatures become larger as a 
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FIG. 4: Spatial variations of the local electron density m for 
various temperatures (a) T=0.0 (b) T=0.11 and (c) T=0.17. 
The strength of the on-site repulsion U = 2.4. The averaged 
electron density is fixed at n = 0.88. 

result of stronger AF order near the vortex core in under- 
doped samples. For U = 2.4, the Q v /e versus tempera- 
ture curve (see the solid line in Fig. 5) exhibits a discon- 
tinuity at T s ~ 0.1 which reflects the characteristic tem- 
perature between two different symmetries. For U = 2.2, 
the spatial distribution patterns for SDW and CDW are 
always isotropic and two-dimensional with Taf ~ 0.09 
where a appreciable discontinuity in Q v /e (see the dot- 
ted line in Fig. 5) is present. 

Recently STM experiments [H HJ showed the exis- 
tence of inhomogeneities in the electron density even for 
the optimally doped E^S^CaC^Og+x, there the hole 
density varies from 0.1 to 0.2. From our previous pa- 



per (2^j and the present study, we found that there may 
exist two types of vortex cores either negatively charged 
or positively charged. By increasing T, the number of 
negatively charged vortices decreases, while positively 
charged vortices increase. This may have profound effect 
on the mixed state transport properties in HTS. Finally 
we wish to point out that our numerical study has been 
performed in a strong magnetic field. So far we are not 
able to extend this calculation to a much weaker field. 
However, according to the extrapolation scenario [2nj | . 
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FIG. 5: The vortex charge number Q v /e for U = 2.4(solid 
line), U = 2. 2 (dotted line) as a function of temperature T. 
The averaged electron density is fixed at n = 0.88. 



the magnitude of the induced AF order follows roughly a 
linear relationship with B. The strength of the AF order 
and the magnitude of vortex charge are estimated to be 
much smaller as B is reduced to the experimental value 
(~ 9 AT). Despite this difficulty, our numerical approach 
should give a consistent description to the behavior of 
vortex charges in HTS. 
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